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Abstract 

We present our software to compute gamma-ray emission from inverse-Compton scattering by cosmic-ray leptons 
in the heliosphere, as well as in the photospheres of stars. It includes a formulation of modulation in the heliosphere, 
but can be used for any user-defined modulation model. Profiles and spectra are output to FITS files in a variety of 
forms for convenient use. Also included are general-purpose inverse-Compton routines with other features like energy 
loss rates and emissivity for any user-defined target photon and lepton spectra. The software is publicly available and 
it is under continuing development. 
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1. Introduction 

The importance of gamma-ray emission by inverse- 
Compton (IC) scattering of cosmic-ray electrons and 
positrons on the photon field of individuals stars was 
first realized in 2006 ||5]|6). The same mechanism for 
emission coming from the Sun was found to be even 
more interesting. In fact, the Sun was predicted to be 
an extended source of gamma-ray emission, produced 
by IC scattering of cosmic -ray electrons on solar pho- 
tons in Q, flU and @. Following the discovery of solar 
emission from the quiet sun in EGRET data [7|, now 
Fermi-LAT is so sensitive that even such weak emis- 
sion can be detected with high significance and studied 
in detail |Q]. Hence, propagation of leptons in the in- 
ner heliosphere can be investigated, which is otherwise 
impossible. Inverse-Compton emission has a broad dis- 
tribution on the sky with maximum intensity in the di- 
rection of the Sun, but extending over the entire sky at a 
low level. Hence, solar IC is also important as a back- 
ground over the sky to be accounted for in studies of 
Galactic and extragalactic gamma-ray emission. 

Calculations of the Galactic IC emission have usually 
assumed a smooth interstellar radiation field, but in fact 



a large part of the Galactic luminosity comes from the 
most luminous stars, which are rare. Therefore we ex- 
pect |8j the interstellar radiation field, and hence the 
inverse Compton emission, to be clumpy. Also the pre- 
dicted IC emission from Cygnus OB2 J6) can be mod- 
eled for example as in (8), and contributes to the the 
emission from the Cygnus region observed by Fermi- 
LAT 0. 

It is important to have a reliable reference IC model 
available for general use. We present here our C++ soft- 
ware to compute IC scattering from the heliosphere, as 
well as the photospheres of stars. It includes a formu- 
lation of modulation in the heliosphere, but can be used 
for any user-defined modulation model. It outputs an- 
gular profiles and spectra to FITS files in table format 
in a variety of forms for convenient use. The software 
is publicly available^ and is under continuing develop- 
ment. 



'http://sourceforge.net/projects/stellarics 
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2. The software package 

The software computes inverse-Compton scattering 
from the heliosphere and the photosphere of individ- 
ual stars. The software uses general-purpose inverse- 
Compton routines with other features like energy loss 
rates and emissivity for any user-defined target photon 
and lepton spectra. It is written in C++ and is modular. 
There are C++ classes for IC cross sections, emissiv- 
ity spectra, stellar radiation fields, electron and positron 
spectra. It can compute both isotropic and anisotropic 
IC scattering. A driver routine is provided and can be 
adapted by the user as required. New models of elec- 
trons, positrons and modulation can be easily added. It 
has an optimized emissivity computation for spectral in- 
tegrations and it can be run in parallel; the computation 
parameters allow a user-defined compromise between 
high resolution and reasonable computation time. It is 
very flexible since it contains user-defined parameters 
(energy ranges, integration steps etc.). Output is pro- 
vided as FITS files (in table format), in various forms: 
angular profiles, spectra and differential and integrated 
flux. Output as idl commands is also provided. It has 
been tested with GNU and Intel compilers, and option- 
ally uses OpenMP for faster execution on multiproces- 
sor machines. A sample output dataset is provided with 
the package to illustrate the format and check the instal- 
lation. 

2.1. Main program 

The software package includes an example 
{solarJcudriver.ee), which contains the parame- 
ters that can be edited by users depending on their 
requirements. This parameters to be defined are: energy 
range and grid factor for both cosmic -ray electrons and 
gamma rays, grid spacing (linear or logarithmic) of 
the angle from the Sun or the star, integration range 
and steps, parameters of the star (temperature, radius, 
distance), electron spectra models including cases with 
free parametrizations. 

2.2. Calculation 

The calculation of the IC emission integrated 
along the line of sight is contained in the class 
SolarlC. It uses the input parameters from 
solar Jc Ariver.cc, the electron spectrum parametriza- 
tion from LeptonS pectrum, the photon field of the star 
or Sun from S tarPhotonF ield, as a black-body source, 
and the IC cross-sections based on the Klein-Nishina 
formula from Inver seCompton. The option to use a 
logarithmic angular grid is useful to resolve the rapid 
variation near the solar disk while still covering the full 
angular range required. 



2.3. Electron spectrum 

The electron and positron spectra can be de- 
fined by the user by adding named models to the 
LeptonS pectrum class. The model parameters are 
freely definable. Several sample models are provided, 
including those based on existing publications OHIO"). 
They will be regularly updated with new models. Mod- 
ulation is described with the force-field approximation, 
but this can be freely extended to more realistic models. 

2.4. Inverse Compton differential cross-section 

Cross-sections are calculated in class 
Inver seCompton. Both isotropic and anisotropic 
Klein-Nishina cross-sections are included, as described 
in |[3j [7]| . The Thompson cross-section is also included 
for a cross-check. The IC emissivity spectrum for a 
given lepton spectrum and radiation field is computed 
here, and this can be used for other applications such as 
interstellar emission. An energy-loss computation for a 
given radiation field is also available here. 

3. Output 

The output is a FITS file with multiple extensions for 
convenient use e.g. directly with HEASARC's fv tool, 
or other plotting tools. Extension names and contents 
of the output FITS file and respectively units are listed 
in Table [Tj Alternatively the FITS file can be read by 
any program using cfitsio, or output as ASCII using 
fv. To plot a differential solar spectrum, as for example 
in Figure [T] run a solar model, choose extension 6 and 
plot the appropriate column against angle. To plot solar 
profiles of integral intensity choose extension 2, and plot 
the appropriate column against angle, as in Figure [2] 

In addition, the program outputs an idl program 
which generates spectra and profiles directly: this op- 
tion has been used to generate the plots in this article. 

4. Examples 

Figs [T] and [2] show the solar spectrum at various an- 
gles from the sun, and intensity profiles above various 
energies, for a typical CR lepton spectrum. For the pro- 
files a logarithmic angular grid was used to resolve the 
behaviour near the solar disk. 

This plot illustrates the extension of the emission to 
the entire sky, and also demonstrates an interesting fea- 
ture: due to the anisotropy of the solar radiation field the 
emission is almost zero towards the solar disk, the emis- 
sion beyond the sun being occulted. This effect should 
be visible in future high-quality Fermi-LAT data, and 
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Table 1 : Extensions names and contents of the output FITS file and corresponding units. I(E, 0) is the differential intensity for gamma-ray energy 
E at angle 6 from the solar centre. 



Extension's Number and Name 


Contents 


Row 


Column 


Units 








Variable 


Variable 






1 . Differential intensity .profile 


I(E,0) 


E 





cirT 2 sr~ 


1 s- 1 MeV 


2. Energy JntegratecLprofile 


/(> E, 0) 


E 





cirT 2 sr~ 


1 s- 1 


3. Angle JntegratecLprofile 


I(E, < 0) 


E 





-2 -1 

cm s 


MeV" 1 


4. Energy _and_angleJntegrated_profile 


/(> E, < 0) 


E 





-2 -1 

cm s 




5. Spectrum_f or .angles 


KE,6) 





E 


cirT 2 sr~ 


1 s- 1 MeV 


6. Spectrum_times_Esquared_for .angles 


I(E, 0) x E 2 





E 


crrT 2 sr~ 


1 s- 1 MeV 


7. Spectrum _for_integrated_angles 


I(E, < 0) 





E 


-2 -1 

cm s 


MeV~' 


8. Spectrum_times_Esquared_for_integrated_angles 


HE, <6)xE 2 





E 


-2 -1 

cm s 


MeV 


9. Energies 


Energies 


E 




MeV 




10. Angles 


Angles 







degrees 






10 2 10 5 10" 
energy, MeV 



10 B 



Figure 1 : Example of output: Intensity spectrum of solar IC at (from 
top) 0.3°, 5", 10", 20", 30°, 40" from the Sun's center. 



can help to distinguish the surface emission (cosmic-ray 
hadronic interactions) from the IC component. 

5. Stellar spectra 

Spectra for stars have been presented in J5]|6l, and 
can be obtained with the present package by defining 
the relevant stellar parameters: luminosity, temperature 
and distance. 




1 10 
angle from sun, degrees 



100 



Figure 2: Example of output: Angular profile of the IC emission as a 
function of the angular distance from the Sun, (from top) > 10 MeV, 
> 100 MeV, > 1 GeV and > 10 GeV. 



6. Conclusions 



We have presented our software for calculating the 
IC emission from single stars and the Sun. The software 
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is under continuing development, taking into account 
updated observations in gamma rays and cosmic rays. 
It is used as a standard model for generating the solar 
IC emission as a background component and it will 
be especially useful for evaluating future data from 
Fermi-LAT. 
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